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This  report  1*  the  result  of  research  performed  in  the  Computational 
Aerodynamics  Group,  Aerodynamics  and  Airframe  Branch,  Aeromechanics  Division, 
Flight  Dynamics  Laboratory  from  June  1984  to  July  1985.  This  report  was 
prepared  by  Dr  Joseph  J.  S.  Shang  and  Mr  Stephen  J.  Scherr,  under  Work  Unit 
2307N611,  "Computational  Fluid  Dynamics." 

Two  separated  presentations  to  the  American  Institution  cf  Aeronautics 
and  Astronautics  Scientific  Meeting  were  derived  from  this  effort .  These 
presentations  are  entitled  "Kavler-Stokes  Solution  of  the  Flow  Field  Around 
a  Complete  Aircraft,"  A1AA  85-0J59-CP  and  "Three-Dimensional  Body-Fitted 
Grid  System  for  c  Complete  Aircraft,"  A1AA  86-0428  respectively. 
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1. 1ST  OF  SYMBOLS 


A  Reference  area  for  aerodynamic  force  coefficient, 

C_  Drag  coefficient,  2D/p  l’  2 A 

I)  w  OO  00 

CT  Lift  coefficient,  2L/p  U  2A 

Cp  Specific  heat  at  constant  pressure 

e  Specific  internal  energy,  C^T  +  0.5(u2+v2+v2) 

F,G,H  Flux  vectors  of  conservation  equations 
Stagnation  enthalpy  of  the  freestream 
h  Enthalpy 

1  Identity  matrix 

i,j,k  Unit  vectors  in  Cartesian  frame 

M  Mach  number 

n  Surface  outward  normal 

p  Static  pressure 

Pr  Molecular  Prandtl  number 

Prt  Turbulent  Prandtl  number 

q  Rate  of  heat  transfer,  q  -  -k7t 

r^  Radius  of  the  blunt  nose,  0.4826  cm 

Rey  Reynolds  number 

St  Stanton  number,  equation  12 

t  Time 

T  Temperature 

u,v,w  Velocity  components  In  the  Cartesian  frame 

U  Dependent  variables  U(p,pu,pv,pv,pe) 

x,y,y  Coordinates  of  the  Cartesian  frame 

1  Specific  heat  ratio 

t  Lddy  viscosity  coefficient 


UST  OF  SYMROLS 
(continued) 


$.n,C  Transformed  coordinates 

W  Molecular  viscosity  coefficient 

p  Density 

T  Stress  tensor 

Subscripts 

*  Freestream  condition 

o  Stagnation  condition 

v  Surface  condition 

Superscript 

T  Transpose  of 


I.  INTRODUCTION 


A  continual  and  systenatlc  progression  of  path-finding  research  in 
solving  the  Navier-Stokes  equations  has  paved  the  way  for  the  first  attempt 
to  simulate  numerically  the  flow  field  around  a  complete  aircraft.1  4  In 
the  present  Investigation,  the  hypersonic  research  aircraft  X24C-10D  is 
selected  for  its  rather  complex  blunt  leading  edges,  canopy,  strake,  fin 
and  wing  formations  which  are  typical  of  those  encountered  in  a  modern 
aerospace  vehicle  (Figure  1).  For  the  chosen  aircraft,  a  detailed  experi¬ 
mental  data  base  also  exists,  in  the  form  of  a  flow  field  pitot  pressure 
survey,  surface  pressure,  heat  transfer,  force  and  moment  measurements. 5* 6,7 
This  experimental  data  establishes  a  solid  criterion  for  assessment  of  the 
accuracy  of  the  simulation.  The  validating  comparison  reveals  technological 
emphases  and  research  opportunities  for  our  further  endeavors  in  computa¬ 
tional  aerodynamics. 

In  principle,  a  composite  solution  using  varying  degrees  of  simplifica¬ 
tion  of  the  Reynolds  averaged  Navier-Stokes  equations  for  the  flow  field 
around  a  complete  aircraft  is  possible.  However,  prior  knowledge  or  unique 
insight  into  the  flow  field®’9  structure  is  required  to  implement  this 
zonal  approach.  Additional  evaluations  of  accuracy  for  each  zonal  solution 
and  the  boundary  conditions  for  contiguous  zones  become  necessary.  As  t> 
benchmark  calculation,  the  uncertainty  is  completely  removed  by  seeking  << 
solution  of  the  Reynolds  averaged  Navier-Stokes  equations.  In  order  to 
depict  adequately  the  Investigated  aircraft  configuration,  a  mesh  system 
consisting  of  475,200  nodes  was  adopted.  At  present,  processing  and 
managing  this  huge  amount  of  data  requires  a  very  large-scale  vector 
processor  with  an  external  storage  device.  Therefore,  a  major  portion  of 
the  present  study  is  concentrated  on  the  adaptation  of  our  computational 


I 

procedure  to  the  CRAY  XMP-22  vector  processor  with  a  16  million  word  MOS  j 

i 

memory  unit,  the  solid  state  storage  device  (SSD) .  In  essence,  the  adapts-  1 

tion  of  a  numerical  procedure  to  a  very  large-scale  computer  comprises 
construction  of  a  data  base  according  to  the  computer  architecture,  and 

i 

implementation  of  the  procedure  so  as  to  minimize  data  movement.1’3  This 

I 

is  the  first  consideration  so  that  one  can  extract  the  maximum  data  pro¬ 
cessing  efficiency  from  a  vector  processor. 

In  order  to  achieve  the  stated  objective,  a  careful  matching  of 
grid-point  system,  numerical  algorithm  and  computer  architecture  is  essen¬ 
tial.  Since  there  are  no  sharp  leading  edges  in  the  Investigated  configura¬ 
tion,  a  composite  grid  system  specially  designed  for  high  local  numerical 
resolution  is  unnecessary.1*  A  single  surface  clustered  body-conformal  grid 
system  is  sufficient  to  facilitate  the  aerodynamic  force  calculations.  The 
adopted  MacCormack  explicit,  unsplit  algorithm10  is  known  to  advance  in 
time  with  the  minimum  number  of  data  accesses  and  has  attained  extremely 
high  vector  efficiency  by  continual  improvement.1’3  Most  importantly,  the 
stencil  of  this  algorithm  requires  no  more  than  four  contiguous  nodes  to  be 
processed  simultaneously  in  any  given  coordinate.  This  feature  leads 
naturally  to  a  page  structure  in  grid  point  system  construction,  anc 
eliminates  the  need  for  interface  boundary  conditions  or  repetitious 
regrouping  of  data  to  accomplish  a  particular  orientation  for  a  numerical 
sweep.11  Finally,  the  minimum  amount  of  data  flow  between  the  external 
memory  unit  SSD,  and  the  central  processor  is  obtained  by  partitioning  the 
data  into  separated  blocks  sized  according  to  the  maximum  number  of  aata 
pages  which  can  be  accommodated  in  core  memory.  These  data  blocks  will  be 
cyclically  rotated  into  the  central  processing  unit  as  required  by  the 
numerical  sweeps  for  convergence. 
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The  numciical  simulation  vas  per  formed  for  the  flow  field  around  the 
entire  aircraft  at  an  angle  of  attack  of  six  degrees  with  a  nominal  Mach 
number  of  5.95  and  the  characteristic  Reynolds  number  of  16. 4  million/m. 

The  benchmark  numerical  results  will  first  be  evaluated  for  accuracy  in 
practical  applications  by  prediction  of  the  surface  pressure  distributions, 
the  shock-wave  structure,  the  heat  transfer,  the  surface  shear  stress 


pattern,  and  the  Integrated  aerodynamic  force  coefficients  whenever  comparison 
data  are  available.  Since  the  entire  flow  field  structure  is  reconstructed 
by  the  numerical  simulation,  the  pertinent  characteristics  of  the  stream 
around  various  components  of  the  aircraft  are  identified  and  highlighted 
tor  further  investigations.  These  include  the  flows  around  the  canopy,  the 
swept-back  blunt  fins,  and  the  induced  vortex  dominated  stream  in  the  lee 
of  the  aircraft.  Then  the  overall  numerical  procedure  is  assessed  for  its 
relative  efficiency,  practical  application  value,  and  the  areas  of  emphasis 
for  future  aerodynamics  research. 

II.  GOVERN IMG  EQUATIONS 

The  three-dimensional,  compressible,  Reynolds  averaged  Navier-Stokes 
equations  in  chain-rule  conservative  form  are  given  as1 


C  m/1!  12  3H.T 

at  ’  34  '  34' 


;  ,  W3F  3G  3H.T 

(,|x’ny’I12)  (3n  ’  3n  ’  dr? 


+  (c  r  r  )(il  12  M)T  _  o 
Ux’VCzM3t  *  3t  ’  3iT  ° 


(1) 


It  is  easily  identifiable  that  the  dependent  variable,  U  consists  of 
p,pu  and  pe 


U  -  U(p,pu,pv,pw,pe) 


(2) 


I 


i 

i 
• 
I 

The  Ilux  vectors,  F,  G,  and  H  are  simply  the  components  of  vectorial 

quantities  of  the  continuity  and  energy  equations,  as  well  as  the  tensorial 

components  of  the  momentum  equations  contained  within  the  divergence  j 

I 

operator  of  the  conservation  law.1 

The  closure  of  this  system  of  equations  is  obtained  by  including  the 

j 

Baldwin-Lomax  turbulence  model12  for  the  Reynolds  stress  tensor  and  by  j 

< 

assigning  a  turbulent  Prandtl  number  with  a  value  of  0.9  for  the  heat  flux, 

I 

x  -  -  [p  +  2/3(y+e) V‘u]I  +  (y+e)  def  u  (3)  ; 

i 

5  -  -  +  !£)  VT  (A)  | 

pr  Prt 

as  well  as  the  equation  of  state  for  a  perfect  gas  and  Sutherland's  formula 
for  the  molecular  viscosity  coefficient.  The  description  of  this  differen¬ 
tial  system  of  equations  is  completed,  once  the  appropriate  and  well-posed 
initial  and  boundary  conditions  are  specified. 

In  the  present  analysis  we  have  taken  advantage  of  the  no-yaw  flight 
condition,  requiring  only  half  of  the  configuration  he  computed  with 
respect  to  the  plane  of  symmetry.  This  simplification  applied  together 
with  the  Idea  of  analytic  continuation  results  in  the  leading  coordinate 
surface  (£“0)  at  the  blunt  nose  becoming  a  developable  surface  by  reflection. 

In  practice,  the  grid  system  is  constructed  by  overlapping  two  initial  data 
pages  in  the  £  coordinate.  The  only  constraints  to  achieve  this  topological 
requirement  are  that  the  basic  cylindrical  coordinates  must  transform  into 
a  spherical  system  near  the  blunt  nose  region  and  the  number  of  peripheral 
nodes  defining  the  snrface  must  be  an  even  integer.  However,  the  latter 
constraint  Is  imposed  only  for  the  purpose  of  convenience.  In  this  manner, 
the  adopted  grid  system  is  body  conformal  and  simply  connected,  and  the 
entire  computational  domain  can  be  mapped  into  a  unit  cube.  Details  of  the 


A 


construction  of  the  mesh  system,  including  the  topology  definition  and  the 
grid  generation  procedure,  nay  be  found  in  an  accompanying  paper13  and  in 
Appendix  A.  In  the  transformed  space,  the  initial  and  boundary  conditions 
are  summarized  as  follows: 

Initial  Condition 

U(0,C.n,O  -  (5) 

The  5=0  surface  is  purely  a  numerical  boundary  created  by  the  grid  system. 
Since  the  two  initial  data  pages  in  5  are  mirror  images,  the  boundary 
values  will  be  updated  automatically.  The  5=1  surface  represents  the 
downstream  boundary.  The  extrapolation  condition  is  commonly  used 


355=1 


On  the  body  curface,  n=0,  the  no-slip  condition  for  velocity  components 
and  the  surface  temperature  duplicating  the  experimental  condition  are 
given.  The  value  of  density  on  the  solid  contour  is  derived  from  the 
isobaric  condition  normal  to  the  surface. 


u=!v=w=0 


T  =  316. 7°k 
w 


n* Vp  =  0 


where  n  is  the  n=0  surface  outward  normal 


n  1+n  j+n  R 
x  yJ  z 

/  n“-+r i2+n2 
x  y  z 

i,j,k  denote  the  unit  vectors  in  the  fundamental  Cartesian  frame.  The  far 


field  boundary  has  transformed  into  the  n=l  surface.  This  boundary  is 


located  outside  the  enveloping  bow  shock  wave  system,  where  the  flow 
retains  the  unperturbed  freestream  value 

D  -  U  (10) 

Ou 

The  C=0  and  surfaces  lay  across  the  plane  of  symmetry.  The  condition 
of  symmetry  prevails 

au  (11) 

9C 


III.  NUMERICAL  PROCEDURE 

All  calculations  were  performed  on  a  CRAY  XMP-22  vector  processor  with 
a  16  million  word  MOS  memory  unit,  the  solid-state  storage  device  (SSD). 

The  SSD  appears  logically  as  a  disk  drive  with  random-access  storage.  For 
reasons  of  efficiency,  at  present  its  operational  configuration  is  limited 
to  non-multi-tasked  applications  (no  parallel  processing).  Therefore,  only 
one  central  processor  unit  of  nearly  one  million  words  of  external  memory 
was  available.  A  mesh  system  consisting  of  475,200  nodes  was  utilized  to 
depict  the  complete  aircraft  configuration.  Specifically,  one  hundred  ten 
(110)  nodes  were  assigned  to  define  the  cross-section  of  the  aircraft  in 
one  hundred  twenty  (120)  streamwlse  stations.  Thirty-six  (36)  nodes  were 
used  to  accommodate  the  computational  domain  bounded  by  the  solid  contours 
and  far  field.  In  short,  a  (110,  36,  120)  mesh  point  system  was  adopted. 

In  the  present  code,  only  five  dependent  variables  lT(p,pu,pe),  three 
independent  variables  x,y,z  and  the  eddy  viscosity  coefficient  required 
permanent  storage.  The  predictor  level  of  dependent  variables  and  nine 
components  of  the  coordinate  transformation  derivatives  are  generated  as 
needed.  Therefore.,  the  maximum  amount  of  external  data  storage  is  4,?76,800 
words. 

MacCormnek'8  explicit,  unsplit  algorithm  vectorizes  easily  with  e-ple 
opportunities  for  chaining  operations.  This  reliable,  conditionally 


MacCormack's  explicit,  unsplit  algorithm  vectorises  easily  with  ample 
opportunities  for  chaining  operations.  This  reliable,  conditionally 
stable  scheme,  when  Incorporated  with  a  time-vsrp  acceleration  procedure 
has  demonstrated  efficiency  gains  of  a  factor  of  six  or  more  for  solutions 
which  possess  a  steady-state  asymptote. J ’ 14  For  this  explicit  algorithm, 
only  four  contiguous  nodes  in  any  given  direction  are  required  to  be 
processed  simultaneously,  the  fourth  node  of  data  being  included  to  support 
the  fourth-orcer  smoothing  terms.  The  stencil  permits  a  streemwise  page 
structure  of  data  in  contrast  to  that  of  an  implicit  scheme  which  requires 
alternating  structures  of  the  data  array  to  support  differing  orientution 
of  numerical  sweeps.  For  the  aforementioned  stencil,  one  detects  that  only 
the  four  cyclic  predictor  pages  currently  in  use  are  required  to  be  retained. 
The  predictor  and  corrector  sequence  can  be  performed  within  one  sweep  by 
lagging  the  corrector  calculation  by  two  spatial  pages.  The  coordinate 
transformation  derivatives  are  repeatedly  generated  by  a  three-point, 
second  order  differencing  scheme  as  needed  without  allocating  a  permanent 
storage  space.  Those  additional  arithmetic  operations  substantially  reduce 
memory  requirements  but  yield  an  in-core  data  processing  rate  of  3.077x10  5 
seconds  per  iteration  per  node  point. 

In  essence,  computation  of  the  three-dimensional  Reynolds  averaged 
Navier-Stokes  equations  becomes  an  intensive  data  transfer  proceduie 
between  the  centra]  processor  unit  (CPU)  and  the  data  storage  device  (ESD). 
For  the  utilized  computing  system,  the  most  efficient  transfer  rate  for 
intermediate  size  data  sets  Is  around  a  halt-million  words  (670  million 
bits  per  second).  However,  in  the  present  investigation,  one  hundred 
twenty  (120)  streamwise  stations  were  divided  into  twenty  data  blocks  of 
six  pages  each.  The  minimum  bandwidth  which  can  support  the  four  page  data 
base  is  the  six-page  block.  The  additional  two  pages  provide  the  coordinate 


transformation  derivatives  at  the  four-page  data  interfaces.  F.ach  blork  of 
data  was  assigned  a  data  set  of  213,840  words.  Iwo  blocks,  or  twelve  pages 
oi  data,  were  simultaneously  resident  in  the  CPU  during  the  numerical  sweep 
between  the  upstream  and  the  downstream  boundaries  (Fig.  2).  The  rest  of 
the  in-core  memory  space  was  taken  up  hy  the  temporary  dependent  variable 
arrays  and  the  coding  Instruction.  Individual  data  blocks  were  rotated 
into  the  CPU  in  sequence  as  required  by  the  two-page  legging  cycle  of 
predictor-corrector  procedure.  In  essence,  the  computer  rode  serves  as  a 
data  manager  to  control  the  four-page  predictor-corrector  cycle  and  the 
twelve-page  external  storage  rotation.  A  measure  of  the  relative  efficiency 
for  adaption  of  numerical  procedure  to  a  very  large  computer  with  an 
external  memory  device  is  the  degradation  of  the  data  processing  rate  (DPR) 
in  comparison  with  in-core  performance.  For  the  present  investigation,  a 
total  of  4,276,800  words  processed  by  a  880,000  word  CPU,  the  degradation 
is  less  than  23  percent  in  DPR.  The  data  processing  rate  in  conjunction 
with  the  SSD  is  consistently  3.7814x10  b  seconds  per  iteration  per  grid 
point. 

As  a  pleasant  surprise,  the  page  and  block  data  structure  has  been 
demonstrated  to  reduce  the  required  computing  time  for  overall  convergence. 
Due  to  the  numerical  sweep  bias,  the  upstream  numerical  simulation  tends  to 
converge  sooner  than  that  of  the  downstream.  As  the  data  converged  In  an 
upstream  block  need  not  be  rotated  into  the  CPU,  the  totci  number  or  data 
blocks  can  be  truncated  accordingly.  A  net  saving  in  computing  time  of  up 
to  thirty  percent  has  been  realized  by  suppressing  unnecessary  calculation 
in  the  converged  region.  Computation  performed  ioj  ibis  truncated  dona u 
is  effectively  equivalent  to  a  zonal  method,  but  has  no  need  of  explicit 
zonal  boundary  conditions. 


For  the  numerical  simulation  of  a  complete  aircraft,  the  surface 
normal  grid  spacing  required  to  resolve  the  shear  aliens  or  heat  ttansler 
can  no  longer  he  gouged  by  a  single  selected  value  of  the  law  of  the  veil 
variable,  y*.  Consideration  Bust  be  given  tc  the  large  upstream  ratio 
between  normal  and  peripheral  nesh  sise  and  the  growth  rote  of  the  boundary 
layer  along  the  entire  configuration.  From  two  different  estimates  of  the 
laminar  sublayer  thickness  based  upon  •»  turbulent  boundary  layer  scaling 
law,  the  mininun  surface  normal  step  size  was  determined  to  be  9.1x)(<  1  cm. 
This  nominal  value  was  used  as  a  guide  tc  wet  an  appropriate  oveiw! 1  grid 
spacing  and  s  finer  mesh  step  In  the  nose  region,  where  the  surface  normal 
step  size  is  s  factor  4.09  time  finer  than  the  nominal  value.  The  periph¬ 
eral  nude  distribution  is  weighted  by  the  local  radlut  ot  curvature.  The 
resulting  stretching  is  relatively  limited  in  magnitude,  gemaily  com  it ed 
within  a  fraction  i>f  a  basic  unit.  However,  .'r:  the  streaawlse  dim  lion, 
the  ratio  between  the  iwir imum  and  the  minimum  step  sire  is  ae  high  as 
19.47.  The  high  grid  del  city  zones  were  concentrated  arourd  the  nom 
,eg1  nn  and  the  locatlonc  where  the  canopy,  strakc  and  fin  protrude  into  tl* 
encoming  stream.  The  hundred  twenty  pages  arc  nearly  equally  dlv1dr< 
between  the  forebody  and  the  afterbody.  Analysis  c‘  the  flow  field  lmmedi 
stely  downstream  of  the  blunt  nose  through  the  emergence  of  the  Mich*  bv  a 
parabolized  Na  v  le  r-Stokes  scheme*^  indicated  t  tie  t  the  mavinus.  fcrld  spacing 
rould  be  used  in  this  region.  The  assigned  value  was  i.a4e  cm.  finally, 
the  steady  asymptote  of  the  nnsri  ic«]  solulior  Is  inniMti  *d  to  be  reached 
when  the  1  ,  norm  between  two  consecutive  noit. ’i/ed  pressure  value: 
than  10  .  This  crlterim  n>  much  less  stringent  t  I  an  the  evaluation 

hea  t  transfer  and  is  < onaldertd  suitable  for  our  present  purpeu* 
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IV.  niSClJSSIUW  OF  NUMERICAI  IhSULTS 
The  discussion  of  results  comprises  two  stjoi  portions.  The  first 
group  of  presentations  le  accompanied  by  experimental  measurements  of  nhock 
structure,  and  heat  transfer  data.b>t>  Then  a  specific  comparison  with  a 
numerical  solution  of  the  patahollzed  Navler-Stokes  conations  for  the 
furchody  Is  presented  and  discussed.  Dhcervations  on  the  characteristics 
of  flow  observed  from  the  numerical  simulation  will  alao  be  de! ineared  to 
Identify  ureas  of  future  research  emphasis. 

Fines  most  cf  the  experimental  measurements  were  collected  or  several 
predetermined  streamvlse  stations,  the  format  of  comparison  Is  constructed 
accordingly.  lypiiitl  results  are  given  by  ahsclsees  In  the  form  cf  the 
normalized  art  length.  Thle  length  Is  insiputed  f roei  the  leeward  meridional 
plane  toward  the  windward  t onnterpart  and  scalec  by  the  total  arc  length  o' 
each  Individual  cron  taction.  In  Figure  3  am.  figure  4  a  total  of  eleven 
strsurwlse  stations  are  depit  ted.  These  stations  icpiesent  a  sequential 
evolution  f  rttr  tie  blunt  nose  region  orti.  the  traiiirp  section.  The 
physltal  locations  span  a  range  of  IS  to  141.^3  radi;  f  t  oir  the  c. ordinate 

orlglr  which  Is  1.  'f  tadli  upstream  of  the  blunt  nose,  rv  *  .483  cm.  The 

*» 

♦  r  r.  1  1  agreement  betwett  data  and  numcilcal  refills  Is  excellent.  In 
figure  *.  the  at  rear*  is*  t umptcff Ion  anil  t»-  icbscuuint  expansion  ,<iound 
tit  tntipv  due  t'  the  recessed  atul  llatt  cuk  J  out  ot  v.ird  surlatt  are 
obvious.  lr  fat  I  ■  ;  i  ■  !  set.  t  1  on.. !  plant'.  Mu  t  .  pans  lm  I  r  re  w  1  ndw.t  r  d  1 1- 
'•.■ward  domain.-,  stands  ml  .  The  surf  ate  pressuit  <M  si  t  Ibut  Ions  alst 
exhibit  a  relative’  ritti  t.  omp  1 1  salon  at  Mu  Mwer  uirmt  if  the  fr.iji.oldal 
toss  section.  This  !rr  m  phenomet  ■  i  ’  ears  a  siting  reset!  '■!».«•  to  thet  ul 
the  hlunt  k  it  g  edge  t  a  delta  w  '  i  (  lhe  curtate  pr  e  s  ■  t  r  t  u  tat  i  tin.  flops 

of  the  a  f  t  e  Miocl'v  att  ,  *  .  en  Sv  Mntt  ■»  .  tail  |  •  ,  ssnrf  «;<!.•  t  epr  est  i  t  -  tbe 


result  of  protruding  fins  oi  stroke  into  the  onenmlng  stream.  The  presM.i*- 
distribution  reveals  the  emergence  of  the  csntiel  fin  near  the  origin  oi 
abscissas,  the  middle  tin,  and  the  stake  from  the  afterbody.  Finally,  the 


stake,  due  to  recession,  ceases  to  be  an  exposed  blunt  edge  (x/r  -  141.25). 
Unfortunately,  the  x/r  *  125.75  station  has  the  only  tuo  pres  i n  laps 
ins Lulled  on  the  leading  edge  of  the  stroke.  The  accordance  between  the 
data  and  numerical  issults  is  excellent. 

Figure  5,  three  pressure  distributions  around  the  central  tin  «> 
different  elevations  (z/r  *  20,  25,  30)  are  given.  These  elevations 
correspond  to  distance*  of  11.92  ca,  3.73  re ,  and  5.74  cn  above  t  be  upper 
nodal  surface.1’  The  collapsing  of  all  pressure  distributions  into  a  narrow 
band  reveala  that  It.  spit*  of  all  upstrsan  diet  urbane  as ,  the  leaverd  flow 
field  Immediately  surrounding  tha  central  fin  remains  nearly  uniform  at 
this  noderate  angle  of  attack.  The  naxiaua  deviation  between  data  at.d 
calculations  la  Halted  tc  a  tsv  parcent. 

The  flow  field  survey  data  by  pit  it  pressure  probes*'  is  presented 
together  with  the  computed  impact  pressure  contour  in  Figure  6.  It  is 
clearly  exhibited  that  the  numerical  simulation  duplicates  the  measurements. 
The  definition  of  the  shock  wave  system  sc  captured  by  the  numerical  mean 
Is  comparable  to  the  data  scattering  band . *  As  we  have  mentioned  earlier, 
the  enveloping  shut  V  wave  system  has  attenuated  substantially  li:  the 
leevcrd  region  to  be  detected  as  a  Mach  wave. 

In  Figure  7,  the  heat  transfer  In  form  cf  Stanton  number  distributions 
In  five  st teamwise  stations  (x/r  *108,  117,  IPS,  75,  1H.5G,  141.25)  are 
depicted.  These  data  collecting  stations  are  Identical  to  there  of  the 
surfai e  pleasure  in  the  afterbody  region.  The  Stantoi  number  is  di fined  bv 
Pt  *  S,  '('.“.(H^-h,)  (12) 
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The  numerical  value  of  the  rate  of  heat  tram-fai  was  computad  dlreit’y  lrom 
tha  prolile*  of  spaciflc  internal  energy  lnataad  of  static  temperature  to 
taka  advantage  ot  tha  monotonic  distribution  of  tha  former  ml  j scent  to  tha 
body  surface.  Tha  specific  foiaula  usad  la 
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Tha  coordinate  transformation  derivatives  are  IncluotO  a* 
obtaining  tha  haat  flux  normal  to  lha  body  surface,  q  -  r'k 
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overall  comparison  between  data  and  computational  result*  la  only  fair. 

The  peak  heating  la  consistently  romentrated  around  tha  exposed  blunt 
leading  edge  of  thr  fins  and  atrake.  Since  no  thermocouple  could  he 
installed  at  thane  landing  edges  because  cl  thalr  size,  a  verification  by 
the  analyll<  solution  baeed  upon  lha  reference  temperature  method  for  the 
stagnation  point  was  performed.16  The  analytic  solution  ylaldad  a  Stanton 
number  ol  1.243x10  1  which  1*  within  tan  percent  (10T;  cl  the  numerical 
prediction  locally,  however ,  the  persistent  .<m!  maximum  discrepancy  (.0X1 
batvaan  data  and  the  praaunt  calculations  Is  located  on  the  windward 
Milace  and  over  a  wide  r.u>ge  ot  the  btreamwlse  stations.  '.his  behavior 

let-  been  commonly  detected  in  numerical  simulation*  of  flow  fields  will 
incidence.  Although  t he  local  grid  density  it  teialivelv  high,  t he  til mum 
notmal  gi  id  sparing  dce.s  not  commit  according  to  the  comptissed  botrdary 
layer  thltkt'css.  This  observation  leads  logit  cl ly  to  tin  peculation  that 
further  gild  retirement  Is  probably  more  iiitii.il  that  the  imprt  \  ement  of 
the  turbulence  model.  A  single  ef  fort  ir  this  nv.i  .  uit  area  w.u'd  yield 
some  Invaluable  basic  knowledge.  Ir  all,  the  pi  event  calm 'at  Ions  u  '  ■  w 
the  trend  o?  data  tn  reveeiing  the  peak  heating  <  f  the  exposed  leading 
edges  end  troughs  ot  heel  transfer  Induced  b>  1 1  •  veparated  iiui  uu 


teparat  ion  rones  on  tin-  ’  ceward  surface  anil  i.*oii  region? 
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A  »«r leu  of  rwirbbl*  numerical  solution*.  using  cite  parabolized 
Kavler-Stokes  equations  has  been  performed  foi  the  forebody  of  the  inves¬ 
tigated  complete  configuration  by  Chaussee  et  al.8'9  Their  composite 
numerical  procedure  initiated  a  full  Navler-Stoket  solution  for  the  blunt 
noae  region  which  served  as  the  initial  boundary  condition  of  the  subsequent 
space-war ching  solution  until  the  strake,  then  the  solution  was  terminated. 
Their  experience  suggested  a  very  coaree  at rennets*  step  size  was  permissible 
between  the  canopy  and  tha  baginning  of  the  strake.  Therefore,  only  twenty 
one  (21)  streamwise  stations  and  the  coarsest  mesh  spacing  (A*  r  1.45  cm) 
were  impleawnted  within  this  region  for  the  present  investigation. 

Surface  pressure  and  the  Stanton  number  distribution  comparisons  with 
the  parabolized  Navlrr-Stokes  solutions  at  their  furthest  downstream 
station  (x/r  -  108)  ere  presented  in  Figures  8  end  9  respect Jve.y .  The 
agreement  In  surface  pressure  distribution  is  very  go od .  However,  the 
presence  of  upstream  pressure  propagation  due  to  tha  strake  lr  detectable 
in  the  present  compulation.  The  small  difference  between  numerical  solution* 
is  therefore  attributable  to  the  omitted  interaction  In  PNS  solution  for 
the  protrusion  of  the  strake.  The  difference  it*  the  heat  transfer  datu  is 
more  noticeable  In  the  'reward  surface  beneath  the  roll-up  of  the  vortu:. 
l*if««r  accordant  *  Is  reached  on  the  windward  sutlore  and  both  numerical 
values  under-pred  i  c  t  the  experimental  rui.tureecnt  .  The  strong  dependence 
or  grid  system*,  adopted  in  the  PUS  colt  ulat  ions  was  reported.0  In  order  to 
avoid  confusion,  tidy  data  from  the  cylindrical  coordinate  system  generated 
by  a  hyperbolic  procedure  'which  give  the  best  agreement  with  the  experi¬ 
mental  data)  le  Included  here  for  coaqrarisur. 

The  windward  and  leeward  prewaure  distributions  .-long  the  plane  of 
symmetry  are  given  by  Figure  10.  The  outstanding  feature  is  clearly  that 

I  t 


for  a  hypersonic  lifting  configuration  the  surf nee  pressure  varies  in 
Magnitude  over  severs]  orders.  In  the  nose  region,  the  pressure  attaint 
the  sulnus  value  of  41.8  rinu  of  the  freest ream;  the  minimum  value  Is 
only  about  O.OS  tines  of  the  unperturbed  value  in  the  corners  between  fins. 
The  windward  pressure  distribution  indicates  the  rapid  expansion  free  the 
stagnation  point  and  maintains  a  gradual  compression  throughout  the  windward 
surface  until  the  expansion  Initialed  at  the  trailing  bcactail.  The 
leeward  surface  pressure  undergoes  expansion  from  the  nose,  recompresses 
over  the  canopy  followed  by  the  expansion  and  finally  rises  over  the 
central  fin.  The  agreement  between  the  parabolized  Navicr -Stokes  solution 9 
and  the  present  result  is  very  good  for  the  forebody,  excepting  perhaps 
Feme  mild  fairing  which  may  have  caused  the  undcrpredlctlon  of  the  pressure 
value  over  the  canopy  by  the  parabolized  Navler-Stokes  procedure.  Only  the 
comparison  of  present  results  and  data  Is  existing  for  the  afterbody  of  the 
complete  configuration.  The  overall  agreement  serves  to  validate  the 
preecnt  effort  end  confirm  the  choice  of  grid  spacing  in  the  forebody 
region . 

V.  POST  PROCESSING 

In  thin  section,  a  reconstruction  of  the  entire  flow  field  around  the 
X24C-10D  research  aircraft  is  attempted  to  gain  understanding  <>f  the 
salient  features  of  its  structure.  The  overall  surface  pressure  distribution 
and  heat  tiunsfer  around  the  alrcicit  is  given  together  with  the  immediate! y 
adjacent  surface  stream! ine  to  depict  the  surtace  shear  stress  pattern  and 
to  identify  the  possible  separated  flow  regions,  however,  the  graphic 
presentation  of  the  flow  field  in  form  of  Mach  contours  is  limited  to  a  lev 
cross  sections.  Finally,  the  detailed  numerical  procedure  used  in  evaluat¬ 
ing  the  aerodynamic  forces  exerted  on  the  aircraft  iu  included.  This 


14 


ji  •  v 


y.v-'jw': 


Information  Is  probably  the  most  critical  to  the  aircraft  denier,  and  should 
be  provided  in  the  numerical  simulation  effort. 


The  overall  surface  pressure  on  the  X24C-10D  experimental  aircraft  is 
depicted  iu  a  perspective  view  in  Figure  11.  (Shang  and  Scherr  plate  1). 

A a  we  have  mentioned  previously ,  the  value  of  surface  pressure  spans  a 
range  of  several  decades  so  an  artificially  compressed  scale  (cube  root 
value  of  pressure)  was  used  to  span  the  pseudo  color  spectrum.  As  antici¬ 
pated,  the  high  pressure  regions  are  easily  recognizable  around  the  nose, 
canopy,  leading  edges  of  st rakes  and  fins  the  aircraft  components  exposed 
to  the  oncoming  stream.  The  extensive  expansion  above  the  leeward  surruro 
downstream  of  the  canopy  is  also  apparent,  but  the  lowest  pressure  level  is 
confined  to  the  inner  fin  root  region  of  the  middle  fin.  The  evidence  of 
Interference  on  the  middle  fin,  originating  from  the  strake,  is  reflected 
by  the  intermittent  high  and  low  surface  pressure  value  ever  the  backward 
swept  fin.  This  phenomenon  Is  verified  by  the  experimental  data.  One  of 
the  interesting  revelations  is  that  the  pressure  along  the  side  surface  of 
the  fuselage  is  nearly  constant  until  the  afterbody  contracts  to  form  a 
recessed  boottAll  where  the  flow  expands  toward  the  downstream.  The 
recompresaion  along  the  plane  of  symuetry  to  realign  the  orientation  of 
vortical  flow  above  the  leeward  surface  is  also  clearly  indicated. 

The  characteristics  of  heat  transfer  to  the  surface  are  depleted  by 
the  static  temperature  of  the  immediately  adjacent  nodes  away  from  the 
surface  (Figure  12,  Shang  and  Scherr,  plate  2).  Since  the  body  surface  is 
maintained  at  a  constant  temperature,  this  local  tempetuture  effectively 
dictates  the  intensity  of  heat  transfer  to  the  solid  surface.  The  high 
heating  regions  generally  correspond  to  the  exposed  elements  of  the  aircraft 
analogous  to  the  surface  prennure.  However,  secondary  high  heating  zones 
are  also  Indicated  in  the  inviscid-viscous  interacting  regions  surrounding 
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the  canopy,  adjacent  to  the  plane  of  symmetry  and  the  areas  between  strake 
and  fins  as  well  as  the  domain  bounded  by  the  fins.  The  only  unusual 
patterns  observed  are  the  overexpansion  of  upwash  leading  to  a  lower  heat 
transfer  rate  around  the  leeward  outboard  corner  and  continuing  to  the 
leading  section  of  the  strakc  where  the  afterbody  recedes.  All  these 
observations  are  supported  by  the  detailed  calculations  presented  in  the 
previous  section  cf  discussion. 

The  effort  involved  in  a  complete  reconstruction  of  the  entire  flow 
field  in  a  graphic  presentation  has  proved  to  be  prohibitive  at  present. 

The  basic  difficulty  of  two-dimensional  projection  of  a  three-dimensional 
phenomenon  can,  in  principle,  be  overcome  by  a  holographic  stereogram. 
Unfortunatelj',  the  supporting  facilities  in  long-haul  communication  and 
limiled  capacity  of  graphic  work  stations  in  simultaneous  display  of  data 
still  fall  short.  Tn  order  to  exhibit  the  peculiar  features  of  the  Inves¬ 
tigated  flow  field,  four  cross-scctlonal  Mach  number  contours  are  given  in 
Figure  13  (Shang  and  Scherr,  plate  3).  At  the  streamwise  station  x/r  = 

39.5,  the  induced  shock  wave  system  over  the  canopy  and  the  dispersion  of 
the  bow  shock  envelope  over  the  leeward  domain  of  this  section  cf  the 
forebody  are  apparent.  At  the  last  forebody  cross  section,  the  attenuation 
of  the  bow  shock  wave  by  the  expansion  from  the  windward  to  the  leeward 
domain  becomes  very  clear.  The  shock  wave  generated  by  the  canopy  by  now 
is  completely  dispersed.  Similarly,  the  vortical  formation  near  the 
leeward  plane  of  symmetry  is  also  obvious,  together  with  the  induced  higher 
Mach  number  pockets  above  and  outboard  of  the  vortices.  At  the  next 
station  downstream,  x/r  ■  125.63,  the  bow  shock  wave  enveloping  the  strake 
impinges  on  the  middle  fin  and  exhibits  significant  distortion  in  structure. 
This  interference  is  observed  persistently  over  the  entire  afterbody  even 


downstream  of  the  location  where  the  leading  edge  of  the  .stroke  sweeps 
backward  and  parallel  to  the  axis  of  symmetry  of  the  afterbody.  In  the 
fast  cross  sectional  plane  presented,  x/r  ■  145.  The  bow  shock  wave 
originating  from  the  blunt  nose  continues  to  exist  In  the  windward  domain. 
The  enveloping  bow  shock  waves  over  the  fins  and  the  £ brakes  dominate  the 
flow  field.  Even  through  the  leading  edge  of  the  stroke  has  turned  parallel 
to  the  fuselage  axis,  its  associated  bow  shock  wave  still  retains  sufficient 
strength  to  Interact  with  the  bow  shock  waves  over  the  middle  fin.  Unfor¬ 
tunately,  little  detailed  Information  is  discernible  by  the  relatively 
coarse  supporting  local  mesh  sysLcni.  The  bounding  shock  waves  above  the 
strakes  are  observed  to  merge  with  the  bow  shock  wave  generated  by  the 
blunt  nose  with  attenuation  around  the  lower  fuselage  convex  corner.  The 
bow  shock  wave  over  the  central  fin  also  interacts  with  the  vortical 
structure  in  the  fin  root  region.  However,  the  interaction  is  rather 
limited  in  intensity  and  is  confirmed  1>y  the  collapsing  surface  pressure 
distributions  around  the  central  fin  (Figure  5).  It  is  also  Interesting  to 
note  that  the  high  Mach  number  pockets  induced  by  vortex  and  expansion 
still  persist  even  to  this  streamwlse  station. 

Mach  number  contours  in  the  plane  of  symmetry  are  depicted  In  Figure 
14  to  serve  as  an  additional  cross  reference.  The  bow  shock  wave  system 
originating  from  the  blunt  nose  is  clearly  identifiable.  In  the  windward 
domain,  the  shock  strength  is  maintained  over  the  entire  corf Iguratlon. 

The  expansion  triggered  by  the  retracted  lower  boattail  surface  is  also 
revealed.  In  the  leeward  domain,  the  shock  wave  system  induced  by  the 
canopy,  its  subsequent  dispersion  and  the  following  bow  shock  waves  over 
the  central  fin  stand  out.  The  envelope  of  the  bow  shock  vave  system 
projected  on  the  plane  of  symmetry  yields  a  total  wave  inclination  of  22 


degrees  over  the  iorebody.  This  value  is  in  perfect  accordance  with  the 
experimental  observation  by  a  shadowgraph. 

The  surface  streamlines  prepared  by  P.  Buning  of  NASA  Aires  Research 
Center  is  presented  in  Figure  15.  This  approximate  surface  shear  flow 
pattern  is  constructed  by  the  first  node  velocity  components  away  and 
parallel  to  the  body  surface.  The  upwash  of  the  nose  region  due  to  the 
Incidence  and  the  secondary  flow  structure  immediately  upstream  of  the 
canopy  are  clearly  revealed.  The  surface  streamlines  converging  at  the 
onset  of  the  canopy  and  diverging  over  the  canopy  surface  suggest  a  confined 
and  small  separated  flow  region  embedded  in  the  junction  between  canopy  and 
fuselage.  The  low  local  heat  transfer  rate  seems  to  confirm  this  observation. 
However,  the  separation  and  reattachment  over  the  leeward  surface  downstream 
of  the  canopy  is  without  any  ambiguity.  The  surface  streamlines  also 
reveal  tendencies  to  converge  around  the  junctions  of  middle  fin,  strakc 
and  the  afterbody  implying  incipient  separation.  The  relatively  low  heat 
transfer  distribution  is  in  correspondence  with  this  idea,  but  the  topo¬ 
logical  constraints  of  the  converging  fins  preclude  a  conclusive  interpre¬ 
tation.  An  interesting  observation  can  also  be  made  on  the  lack  of  a 
horseshoe  vortex  structure  around  the  blunt  fins.17  This  may  be  associated 
with  the  swept  back  configuration  of  fins.  The  drastic  difference  in  flow 
field  structure  between  the  vertical  blunt  fin  and  the  swept  back  fins  may 
warrant  a  detailed  and  separate  investigation. 

The  aerodynamic  force  exerted  on  the  complete  configuration  can  be 
easily  derived  from  the  Cauchy's  theorem 


where  the  stress  tensor  is  given  by 
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The  specific  expression  of  the  stress  tensor  has  been  shown  in  Equation 
(3)  which  contains  the  static  pressure  as  a  component  of  the  normal  stress. 
The  net  aerodynamic  force  components  at  a  point  on  the  control  surface,  n  = 
constant,  in  the  Cartesian  frame  are 
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In  the  present  investigation,  the  control  surface  coincides  with  the  body 
surface.  Since  the  present  analysis  omitted  the  base  configuration,  the 
trailing  surface  was  not  included  in  the  computation.  For  a  hypersonic 
vehicle,  the  drag  contribution  from  the  blunt  nose  and  wave  system  should 
overwhelm  that  of  the  base  regime.  The  present  simplification  will  serve 
as  a  verification  of  this  conjecture  for  engineering  applications.  However, 
the  present  formulation  is  derived  to  be  independent  of  this  approximation 
and  to  retain  the  general  validity. 

The  surface  integral  of  an  arbitrary  plane  defined  by  x(C,0»  y(£,C) 
and  z(£,C)  islfc 


//Fi(x,y,z)  [(x|  +  y|  +  z2)(x2  +  +  z*) 

-  +  y?yc  +  z?Z(;)2]i/2dWC  (16) 
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The  above  expression  can  l>e  further  simplified  by  the  definition  of 
the  Jacobian  of  the  coordinate  transformation  and  some  algebraic  manipula¬ 
tion  to  yield 

//F^x.y.zJKnJ  +  T)*  +  n*)/J 2]  1/2d£,dC  (17) 

Finally,  ve  have  the  net  aerodynamic  force  components  exerted  on  the 
complete  configuration 
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where,  the  Jacobian  of  coordinate  transformation  is  given  by 


The  aerodynamic  force  components  along  any  desired  direction  can  be  obtained 
by  a  simple  projection. 

Tn  the  present  formulation,  the  advantage  of  a  body  conformal  coordinate 
system  is  transparent.  One  recognizes  the  necessity  of  a  zonal  construction 
of  the  grid  system  for  complex  configurations,  however,  consideration 
should  also  be  devoted  to  the  ease  oi  post  processing. 

The  specific  comparison  of  the  computed  drag  and  lift  coefficients  to 


the  experimental  data  is  presented  in  Table  1.  The  difference  in  lift 
coefficients  is  merely  A. 72  percent  end  the  drag  coefficient  has  a  discrep¬ 
ancy  of  6.70  percent.  However,  the  compensating  errors  in  numerical 


underprediction  yields  a  L/D  ratio  within  2.16  percent  of  the  experimental 
measurement.  This  agreement  is  confined  within  the  accepted  experimental 
error  band.  Indirectly,  this  comparison  seems  to  Justify  the  present 
omission  of  the  base  drag  contribution. 

VI.  CONCLUSIONS 

The  first  numerical  simulation  of  the  flow  field  around  a  complete 
aircraft,  the  hypersonic  research  aircraft  X24C-10D,  utilizing  the  Feyuolds 
averaged  Navier-Stokes  equations  has  been  accomplished,  yielding  accurate 
prediction  of  detailed  flow  properties.  The  numerical  results  obtained  at 
a  Mach  number  of  5.9i»  also  indicate  an  excellent  agreement  in  aerodynamic 
force  coefficients  with  those  of  experimental  data. 

The  page  data  structure  collected  into  separated  blocks  and  cyclically 
rotating  from  the  external  storage  device  into  the  central  processor  has 
been  demonstrated  as  an  efficient  and  viable  adaptation  of  the  numerical 
procedure  to  a  very  large  class  computers  without  sufficient  dynamic 
memory.  The  cyclic  data  block  structure  has  also  been  found  to  he  a 
convergence  enhancement  device  by  selectively  processing  data  according  to 
its  rate  of  convergence. 

The  parabolized  Navier-Stokes  solution  generated  for  the  forebody  of 
this  aircraft  has  been  verified  by  the  present  result.  A  composite  numerical 
procedure  for  the  entire  configuration  may  be  desirable  tor  efficiency 
improvement.  Urgent  research  emphases  still  exist  in  continuing  improvement 
of  numerical  efficiency,  long-haul  communication  and  graphic  support  of  the 
computational  aerodynamics  discipline. 
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APPENDIX  A 


Construction  of  the  Grid  Syst  — 

The  grid  system  was  constructed  in  distinct  regions  which  were  then 
merged  to  form  one  computational  block.  The  primary  division  was  between 
the  spherical  grid  at  the  vehicle  nose  and  the  cylindrical  structure  over 
the  rest  of  the  aircraft  body.  Each  plane  of  the  cylindrical  Bystew  was 
further  subdivided  into  three  region*:  a  hyperbolic  grid  next  lc  the  body, 
and  two  transflritely  Interpolated  regions  which  continued  the  grid  fror 
the  hyperbolic  part  to  the  prescribed  outer  boundary  (Figure  16).  This 
arrangement  allowen  *1 '  the  conditions  specified  for  the  grid  to  be  satisfied. 
The  multiple-plane  structure  of  the  mesh  system  was  conveniently  implemented 
by  letting  a  £  ■  constant  plane  alao  be  an  x  «  constant  plane.  The  grid 
was  made  orthogonal  in  each  plane  at  the  vehicle  surface  by  use  of  a 
suitable  hyperbolic  grid  generator.  The  restrictions  on  the  outer  portion 
•>f  the  mesh  were  satlslied  by  the  creation  of  the  composite  structure. 


Smoothness  of  derivatives  between  regions  of  the  grid  was  enhanced  by 
application  of  a  Poisson  operator. 

Hyperbolic  Prid  Generation 

Orthogonality  of  the  grid  at  the  vehicle  surface  was  the  first  require- 
p>eot  for  the  me6h  system.  This  constraint  was  slightly  relaxed  to  require 
only  that  it  be  orthogonal  in  each  x  ■  constant  cross  section.  This  was 
guaranteed  by  use  of  a  hyperbolic  grid  generation  algorithm  originally 
proposed  by  Steger  and  Chausaee19  and  further  modified  by  Kinsey  and 
F>arth.2®  This  algorithm  takes  the  definition  of  grid  orthogonality 


yy  +zr  ^O 

n  t,  n  c. 


(1) 


and  the  transformation  Jacobian 


to  devise  a  b rheme  mapping  vy.s)  to  ( »  .  t ) .  locally  ilnearlrlng  about  the 


point  and  noting  that  an  area  Integrand  c«n  he  written  dV  -  dv  a* 

1/J  dn  dt,  (1)  end  (2)  become 


AJt  ♦  BR  -  F 

C  r 


where  A  * 

ft 

J’, 

c 

£ 

? 

.  B  - 

-i* 

n 

♦yc 

n 

• 

1  - 

0 

V+V° 

.  R 


V  l  ,  .md 


Inspection  of  b  ‘A  shows  that  this  system  ia  hyperbolic  mu]  <  *n  be  marched 
In  n. 


The  finite  difference  solution  scheme  of  the  origins!  algorithm  U 
centrally  cilferenced  in  C  and  backward  differenced  In  r  .  .  t  can  be 

written  as 


♦f  (VLAt)‘kL 
e  K  k  k  ,  j 


where  the  final  term  adds  explicit,  fourth-order  dissipation  as  describee 
In  beam  and  Warning.'"  To  insure  smoothness ,  additional  tempo.*,  and 
special  dlsalpatlon  is  added.  Spaclal  dissipation  corns  from  the  use  ot  a 
more  general  clast,  of  integration  In  r. : 
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When  F“l,  the  original  backward  differenced  scheme  Is  obtained.  lor  d1.'), 
the  error  term  has  a  dissipative  effect  in  the  n  direction.  Revilting  the 
algorithm  in  delta  form  and  adding  tempoial  dissipation  through  an  Implicit 
term,  which  does  not  formally  degrade  the  accura* v  of  the  algorithm,  tlu> 


final  fora  of  th*  hi  Iimh  Is 

[,M|(wkS)^er>A.fi(ij4l-R  ,^»-‘i«v;4lMi-i.)vji-ii(vA)^j 

Not*  that  tliif*  ich«M  require*  th*  l ?-dlm*nslotial)  cell  volumes  at  the  J 

and  1 ♦ 1  levels.  lb*  j  level  volume  ran  be  computed  from  the  *:<  sting  grid 

In  term*  of  the  metric  functions,  using  <?).  Th*  volume  st  i-  estimated 

by  referring  to  sn  average  of  cell  volumes  from  r wo  Ideal  1: to  polar  grids. 

both  these  gride  are  formed  about  circlet  whose  circumference  equals  the 

*r<  1 ei.gth  of  the  body  cross  section  end  both  have  an  eponential  stretching 

in  the  radial  direction.  The  polar  grids  vary  in  their  spacing  In  th. 

<  1 1 > usl erent is  1  direction,  one  being  equl spaced  ana  th*  other  having 

•'pacing  proport  lone'  to  the  original  distribution  of  points  on  the  lodv. 

^b*  j  ♦  i  cell  volume  Is  chose  to  transition  smooth' y  1  roe  the  c orreepond 1 ng 

volume  in  the  unequal!)  spaced  grid,  when  th*  cel'  Is  at  th*  body,  to  the 

olume  In  the  equlsp.u  «c  grid,  when  the  cell  is  ter  from  the  bodv. 

The  petemeters  have  been  used  in  generating  the  grid  control 

dissipation,  spacing,  and  orthogona  1  f t >  .  I  relaxes  the  orf hogora  1  lty  of 

the  grid  as  !  moves  away  f roai  1.0.  This  helps  the  grid  to  adapt  to  the 

- i d*  variation  In  curvature  of  the  body,  without  forcing  it  to  kink  or 

overlap  itself.  Orthogonality  is  forced  at  the  body  surface  by  setting  the 

lollovtng  p«rflih»iufc  at  j-2:  f-1.0,  »  "0.0,  and  t^-U.O.  These  parorciers 

thet:  !ik  ream  lirearly  with  to  their  maximum  value  at  '”1j*  ant*  sf«'>y 

constant  at  that  value.  In  ordei  that  the  radial  llnet.  do  not  spiral  back 

In  toward  the  hoiK  ,  the  value  of  h  Is  returned  to  1.0  between  1*1.  and 

,  .  The  other  values  used  to  genet.itp  the  hvpeihollc  grid.*  were 

IWJI  X 

I'  1.0,  ft  )  “0.0S,  ft.)  -0.04b,  j  -8,  and  ]  «'t0.  The  radial  spacing 

max  e  me x  1  ma >  'I  .  Y  ' 

t  the  wall  was  0.028V  if,  and  the  nominal  total  radiel  distance,  used  in 

estimating  tell  volun.es,  was  2S  r  ♦O.llx-x  ). 

b  c 


The  dfaeribur ion  ot  points  on  the  body  w<>b  controlled  by  use  of  m  t 
length  weighted  by  curvature.  The  hyperbolic  grid  generator  marches  the 
grid  out  from  some  predetermined  body  contour  with  some  given  point 
spacing.  Ar  existing  program  was  used  which,  giver  an  axial  location  x, 
would  produce  a  sequence  ot  points  describing  the  body  cross  section  ,.t 
that  particular  station.  Ihc  curvature,  war  computed  at  each  point  ana 
the  weighted  arc  length,  S(g)-/  l+>nax(r ,0)  ds  was  calculated  along  the 

K 

MX 

curve,  where  1  ie  a  weighting  tactor  for  the  curvature.  Points  were  then 
placed  on  the  curve  no  as  to  be  equlspeced  in  S.  This  function  tenac  to 
rlnetei  points  at  convex  cornets  end  give  a  uniform  distribution  along 
straight  lines  end  at  concave  corners.  This  kind  of  clustering  In  that 
tvpliflly  required  in  an  external  problem.  It  was  also  airanged  that  both 
end  pairs  of  poinit  would  straddle  the  centei 1 lne  of  the  body,  foi  use  In 
Implementing  the  *vwmetr)  «c .tidition.  The  parameter  -*  was  set  to  7er o  for  b 
>  <5.  This  terred  the  points  to  be  equispared  In  oil  length  on  the 

forebody.  a  Increased  linearly  with  c  to  a  value  of  2 .  S  sf  i  *■  >U  and 
lemalned  at  that  level  tor  1  70.  It  wan  K  ped  that  this  method  o' 

piai  4  ny  the  points  would  lest1  to  smooth  changes  between  planes  of  the  griu. 
Composite  Grid  St  rue ture 

The  second  requirement  on  the  »e  ;b  system  was  that  it  confott  u>  a 
smo>  t I  outer  boundary.  An  attempt  was  ir,  ado  t  o  use  vai'lng  cell  volumes  in 
the  hyperboli«  pile,  these  cells  volumes  being  determined  bv  Itet.ting  the 
grlu  generator.  Tt  l  s  attempt  was  unsuccess  t  u  1  ,  an  it  ;  ed  to  .i  kinke<! 
boundary,  and  the  implementation  of  the  r>cthod  proved  to  make  i  •  ••.  tiul  ot 
the  bound  j  i  •-  t  <  o  d  1  f  f  1 i  ul  t  .  Thus  ,  i  t  a  1 t ernat •  method  wa:  t  r  p  1  oved  .  i  *  l  nq 
transflulti  interpolation  from  »  npccitieu  level  ot  the  byperboln  ^tld 
through  p  r  intermeo  ,.'t  <■  smooth  boundary  anc  c>n  to  the  denoted  outei  i.ondarv. 
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The  analytic  description  of  the  outer  boundary  wnr  formulated  based  on 
experimental  results.  A  pitot  pressure  survey  at  the  last  forcbody  station 
(x/rv  ■  128.0)  indicated  the  presence  of  a  bow  shock  wave  of  approximately 
bi-elliptical  shape  (Figure  6).  A  shadowgraph  picture  taken  of  the  experi¬ 
mental  model  showed  that  the  shock  wave  projection  in  windward  domain  grew 
linearly  with  x  at  an  angle  of  0  ■  13.5°.  Thus,  the  outer  boundary  was 
defined  to  he  a  bi-elliptical  shape  centered  at  (0,0)  with  y-axls  length 
A(x)/r^  -  [ 32. 5*(x-128.0  r^)  tan©]*,  upper  z-axls  length  B^x)/!^  * 

[50*(x-128.0  r.)  tan  (G+a)]*,  and  lower  z-axls  length  B„(x)/rk  * 

N  a  n 

[  17.5+(x-128.0  r  )  ten  (0-a)]*.  The  inclusion  of  the  angle  of  attack,  u, 

N 

enables  the  grid  generation  program  to  have  a  wider  range  of  applicability. 

The  multiplicative  factor  *»1.25  increases  the  size  of  the  grin  to  make 

certain  that  the  entire  bow  shock  is  captured. 

In  order  to  match  the  cylindrical  grid  generated  around  the  body  with 

the  apheT i ca 1  grid  generated  around  the  nose,  the  definition  of  the  outer 

boundary  was  altered  at  the  extreme  forward  part  of  the  vehicle.  The 

center  of  the  hemispherical  nose  is  located  at  (x  ,y  ,z  ).  The  line 

e  t  c 

determined  by  B, (x)  ,  the  upper  edge  of  the  grid  system,  was  projected 
forward  all  the  way  to  x*x  .  Thus,  the  radius  of  the  spherical  grid  was  r ^ 

■  b.  (x  ) -z  .  Between  x=x  and  x*x  ♦  20.  75  r..,  an  arbitrary  point  chosen 

iff  c  c  N 

orward  of  rhe  canopy,  all  the  parameters  of  the  bl-elliptical  structure 

varied  1  Inear  1 y — the  center  from  (y  ,7  )  to  (0,0),  and  the  axis  lengths 

c  c 

between  r  and  their  values  at  x  *20. 75  r  .  This  arrangement  produced  a 
t>  c  N 

.mooth  variation  from  the  bl-e  1 1 ipt 1 rn 1  structure  to  a  circular  structure 

at  > * x  . 

c 

An  intermediate  boundary  between  the  outer  boundary  and  the  edge  of 
the  hyperbolic  mesh  was  prescribed  In  order  to  ensure  smoothness  of  the 
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grid  in  the  region  of  the  bow  shock  wave.  This  boundary  was  chosen  to  have 

a  bi-el] iptical  shape,  so  that  all  the  circumferential  lines  between  it  and 

the  outer  boundary  would  be  similarly  shaped.  Thus,  all  the  transition  ! 

between  the  uneven  shape  of  the  hyperbolic  grid  and  the  bi-elliptical  shape 

i 

of  the  outer  grid  took  place  in  the  central  region  (Figure  16) ,  where  the  j 

need  for  numerical  smoothness  was  not  as  pronounced.  The  intermediate 

boundary  was  positioned  to  be  a  fraction  p  *  . 30  of  the  distance  between 

the  outer  boundary  and  a  shape  enclosing  the  hyperbolic  grid.  The  boundary 

of  the  hyperbolic  grid  was  interrogated  to  determine  the  points  of  greatest 

e:;tension  in  y,  (y^.z^),  highest  extension  in  z,  (y2»z2)  and  the  lowest 

extension  in  z,  (y^,z^).  The  center  of  the  bi-ellipse  was  placed  at  (0,Zj) 

with  an  upper  z-axis  length  Bj(x)  *=  z^(l-p)  +  Bj(x)p-Zj,  and  a  lower  z-axis 

length  B^fx)  «  -z^(l-p)  +  B2(x)p  +  Zj.  The  y-axis  length  was  A' (x)  = 

Vj(l-p)  ♦  A(x) p  for  x  >  xc  +  20.75  r^  and  was  A'(x)q  +  l/2(Bj(x)  +  B2(x)(l-q)) 

for  x  <  x  +  20.75  r  ,  where  q  *  c  This  ensures  that  the  inter- 

‘  20.75  rN 

mediate  houndary  approximates  a  circle  afe  x  approaches  x^ . 

The  next  step  in  performing  the  transfinite  construction  was  determina¬ 
tion  of  the  placement  of  points  on  the  intermediate  and  outer  boundaries 
and  specification  of  the  stretching  along  the  radial  lines.  Points  were 
placed  on  the  outer  boundary  so  that  the  arc  lengths  between  them  were 
proportional  to  the  arc  lengths  on  the  edge  of  the  hyperbolic  grid.  Hie 
points  on  the  intermediate  boundary  were  placed  in  a  similar  fashion  with 
an  additional  numerical  averaging.  If  s(n)  is  the  accumulated  arc  length 
through  the  nC^  point,  proportional  to  the  spacing  on  the  hyperbolic  grid 


boundary,  ifi  the  average  arc  length  between  points,  and  o  =  0.1  is  s 

smoothing  parameter,  then  the  actual  arc  length  at  the  point  is  s'(n)  - 

s(n-l)  ♦  [s(n)-6(n-l) ]  ( 1— 0 )  +  s  0.  This  smoothing  relaxes  the  point 

avg 


distribution  so  that  it  approaches  an  equispaced  distribution.  Point 
placement  was  determined  Lot  each  radial  line  separately.  In  the  region 
between  the  hyperbolic  grid  ard  the  Intermediate  boundary,  spacing  was 
determined  by  the  inverse  hyperbolic  sine  iunction  of  Vinokur. The 
spacing  at  the  inner  boundary  coincided  with  the  outer-most  spacing  of  the 
hyperbolic  grid,  and  the  spacing  at  the  intermediate  boundary  was  chosen  to 
be  the  maximum  of  0.125  r^  and  the  spacing  required  for  an  eqtiispaced 
distribution  along  the  radial  line  between  the  inner  and  outer  boundaries. 
Tn  the  other  region  of  the  grid,  an  exponential  stretching  was  used  which 
matched  the  stretching  at  the  intermediate  boundary. 

Generation  of  a  f  '  constant  plane  was  completed  by  application  of  a 
Poisson  averaging  operator  to  the  grid  points  between  j**22  and  j=32,  a 
domain  which  Included  the  entire  middle  region.  This  was  done  to  promote 
smoothness  of  derivatives  across  the  composite  grid.  Some  typical  grids 
can  be  seen  in  Figure  17. 

Spherical  Grid  at  the  Vehicle  Nose 

The  cylindrical  coordinate  system  of  the  mesh  was  easily  transitioned 
into  a  spherical  system  on  the  hemispherical  nose.  Surface  points  were 
spaced  equiangularly  both  latitudlnally  and  longitudinally,  with  radial 
stietching  controlled  by  a  line  from  the  i=7  plane,  which  abuts  the  nose 
region.  As  mentioned  above,  the  parameters  ior  the  composite  grid  are 
varied  so  that  mesh  in  the  i=7  plane  is  nearly  circular.  The  grid  points 
in  the  i*=l  plane  are  the  reflection  in  the  z =z  plane  of  the  points  in  the 
i=2  plane.  Analytic  continuation  of  the  field  variables  in  the  £  direction 
indicates  that  P(k,j,l)  and  P(k,j,3)  are  neighbors  of  P(k,j,2)  spatially 
(Figure  18).  Due  to  the  bilateral  symmetry  of  the  flow,  all  the  scalar 
variables  are  equal  at  P(k,j,l)  and  P <  1 1 1  -  k,j,2),  with  the  velocity 
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vectors  being  equal  in  magnitude  but  reflected  in  tlic  y*0  plane.  Since  the 
spherical  grid  is  equispaced,  the  flow  variables  need  not  be  interpolated 
to  set  boundary  conditions  at  P(k,j,l).  This  type  of  structure  can  be  used 
for  arty  bilaterally  symmetric  flow.  The  singularity  along  the  axis  for  the 
cylindrical-spherical  system  is  physically  avoided  by  the  grid  and  boundary 
tonditions  are  set  across  the  axis  of  symmetry.  One  disadvantage  of  the 
construction  is  that  the  grid  spacing  along  the  surface  of  the  nose  is 
actually  less  than  the  smallest  radial  spacing,  due  to  the  large  number  of 
peripheral  points.  This  is  undesirable  computationally,  but  unavoidable  in 
a  body-fitted  system.  The  extra  work  required  in  processing  was  reduced  in 
the  computation  by  taking  advantage  of  the  data  structure  and  only  computing 
in  the  nose  region  when  necessary. 

In  summary,  construction  of  a  C  =  constant  plane  of  the  mesh  began 
with  the  generation  of  a  hyperbolic  grid  next  to  the  surface.  A  composite 
grid  was  formed  using  transfinite  interpolation  irom  the  edge  of  the 
hyperbolic  grid  through  two  outer  boundaries.  lastly,  a  single  pass  of  .1 
Poisson  smoothing  operator  was  performed.  All  facets  of  the  construct  ion 
were  chosen  to  promote  numerical  smoothness  of  the  grid  and  in  the  ensuring 
computation.  The  mesh  about  the  blunt  nose  v/as  creased  as  a  spherical 
gria,  juxtaposed  to  the  forebody  grid.  Computational  boundary  conditions 
were  provided  for  by  setting  op  the  grid  to  take  advantage  of  analytic 
continual  ion  and  symmetry. 


DATA  STRUCTURE  AND  DATA  FLOW 
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